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ABSTRACT 

We report observations of the flickering variability of the recurrent nova RS Oph at 
quiescence on the basis of simultaneous observations in 5 bands {UBVRI). RS Oph 
has flickering source with [U - B)o = -0.62±0.07, {B-V)o = 0.15±0.10, {V-R)o = 
0.25 ±0.05. 

We find for the flickering source a temperature Tfi w 9500 ±500 K, and luminosity 
Lfi ~ 50 — 150 Lq (using a distance of d = 1.6 kpc). 

We also find that on a (J7 — B) vs {B — V) diagram the flickering of the symbiotic 
stars differs from that of the cataclysmic variables. The possible source of the flickering 
is discussed. 

The data are available upon request from the authors and on the web 
www.astro.bas.bg/- rz/RSOph.UBVRI.2010.MNRAS.tar.gz. 
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1 INTRODUCTION 

In the symbiotic recurrent nova RS Ophiuchi (HD 162214), a 
near-Chandrasekhar-mass white dwarf (WD) accretes mate- 
rial from a red giant companion (e.g., Hachisu & Kato 2001; 
Sokoloski et al. 2006). It experiences nova eruptions approx- 
imately every 20 yr (Evans et al. 2008), with the most recent 
eruption having occurred on 2006 February 12 (Narumi et 
al. 2006). Fekel et al. (2000) found that RS Oph has an or- 
bital period of 455 days and give red giant and white dwarf 
(WD) masses of 2.3Mq and close to 1.4M0 respectively with 
a separation between the components of a = 2.68 x 10^'^ cm. 
For the range of spectral types suggested (Worters et al. 
2007) for the red giant in the RS Oph system, its radius 
is smaller than its Roche lobe, and accretion onto the WD 
may occur only from the red giant wind. 

The flickering (stochastic light variations on timescales 
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of a few minutes with amplitude of a fewxO.l magnitudes) 
is a variability observed in the three main types of binaries 
that contain white dwarfs accreting material from a compan- 
ion mass-donor star: cataclysmic variables (CVs), supersoft 
X-ray binaries, and symbiotic stars (Sokoloski 2003). The 
flickering of RS Oph has been detected by Walker (1977). 
The systematic searches for flickering variability in sym- 
biotic stars and related objects (Dobrzycka et al. 1996a; 
Sokoloski, Bildsten & Ho 2001; Gromadzki et al. 2006) 
have shown that among ~ 200 symbiotic stars known, only 
8 present flickering - RS Oph, T CrB, MWC 560, Z And, 
V2116 Oph, CH Cyg, RT Cru, o Get and V407 Cyg. 

Here we investigate the flickering variability of RS Oph 
in the UBVRI bands and discuss its possible origin. 



2 OBSERVATIONS 

On the night of 2008 July 6, we observed RS Oph simulta- 
neously with four telescopes equipped with CCD cameras. 
The 2m RCC telescope of the National Astronomical Obser- 
vatory Rozhen equipped with a dual channel focal reducer 
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Table 1. CCD observations of RS Oph. In the table are given as follows: the telescope, band, UT-start and UT-end of the run, exposure 
time, number of CCD images obtained, average magnitude in the corresponding band, minimum - maximum magnitudes in each band, 
standard deviation of the mean, typical observational error. 



telescope 


band 


UT 


Exp-time 




average 


min 


-max 


stdev 


err 






start-end 


[sec] 




[mag] 


[mag 


-[mag] 


[mag] 


[mag] 


2008 July 6 




JD2454654 
















2.0 m Rozhen 


U 


19:27 - 21:26 


300 


20 


12.546 


12.293 


- 12.721 


0.103 


0.011 


50/70 cm Schmidt 


B 


19:46 - 21:29 


100,60,20 


50 


12.471 


12.380 


- 12.615 


0.051 


0.005 


2.0 m Rozhen 


V 


18:54 - 21:29 


30 


196 


11.215 


11.079 


- 11.333 


0.052 


0.007 


60 cm Rozhen 


R 


19:36 - 21:33 


15,20 


301 


10.277 


10.172 


- 10.362 


0.036 


0.002 


60 cm Belogr 


I 


19:42 - 21:22 


60 


74 


8.990 


8.855 


- 9.098 


0.057 


0.100 


2009 July 21 




JD2455034 
















50/70 cm Schmidt 


U 


20:58 - 21:52 


120 


21 


12.020 


11.863 


- 12.166 


0.092 


0.013 


60 cm Rozhen 


B 


20:17 - 21:38 


40 


70 


12.030 


11.876 


- 12.277 


0.098 


0.007 


60 cm Belogr 


V 


20:12 - 21:31 


30 


81 


11.009 


10.891 


- 11.236 


0.092 


0.004 


60 cm Belogr 


R 


20:12 - 21:31 


10 


80 


10.219 


10.121 


- 10.406 


0.074 


0.003 


60 cm Rozhen 


I 


20:17 - 21:38 


10 


70 


9.172 


9.070 


- 9.325 


0.066 


0.003 


2009 July 23 




JD2455036 
















2.0 m Rozhen 


U 


19:23 - 21:14 


120 


41 


11.659 


11.482 


- 11.845 


0.081 


0.028 


50/70 cm Schmidt 


B 


20:53 - 21:21 


30 


46 


11.978 


11.857 


- 12.099 


0.062 


0.007 


2.0 m Rozhen 


V 


19:24 - 21:15 


10 


467 


10.932 


10.766 


- 11.093 


0.070 


0.008 


60 cm Rozhen 


I 


19:40 - 21:17 


5 


766 


9.096 


8.970 


- 9.216 


0.050 


0.005 


60 cm Belogr 


I 


19:55 - 21:24 


10 


404 


9.118 


8.998 


- 9.206 


0.048 


0.003 



2009 July 21 " 
I 




R 



B 



I . , \ , , \ , , \ . 1 12.3 I \ , , \ , , \ , , L 
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Figure 1. Variability of RS Oph in the UBVRI bands on 2008 Figure 2. Variability of RS Oph in the UBVRI bands on 2009 

July 6. July 21. 
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Figure 3. Variability of RS Oph in the UBVI bands on 2009 
July 23. 



observed in U and V bands. In the U band a Photomet- 
ries CCD (1024x1024 px, field of view 4.9'x4.9') has been 
used, and in the V band a VersArray 1330B (1340x1300, 
6.3'x6.3'). The 60 cm Rozhen telescope observed in the R 
band (equipped with a FLI PL09000 CCD with 3056 x 
3056 pixels and 5.7'x5.7'); the 50/70 cm Schmidt telescope 
of NAO Rozhen in the B band (SBIG STLllOOOM CCD, 
4008x2672 px, 16'x24'), and the 60 cm telescope of the Bel- 
ogradchick Astronomical Observatory in the / band (SBIG 
ST8 CCD, 1530xl020px, 6.4'x4.2'). 

On the night of 2009 July 21, we observed RS Oph si- 
multaneously with 3 telescopes. The 50/70 cm Schmidt tele- 
scope observed in the U band, the 60 cm Rozhen telescope 
- repeating B and / bands, and the 60 cm Belogradchick 
telescope - repeating V and R bands. 

On the night of 2009 July 23, we observed RS Oph 
simultaneously with 3 telescopes. The 2m RCC telescope 
observed simultaneously in U and V bands, the 50/70 cm 
Schmidt telescope - in B, and the 60 cm Rozhen and Belo- 
gradchick telescopes - in 7 band. 

All the CCD images have been bias subtracted, flat 
fielded, and standard aperture photometry has been per- 
formed. The data reduction and aperture photometry are 
done with IRAF and have been checked with alternative 
software packages. The comparison stars of Henden and Mu- 
nari (2006) have been used. 

The results of our observations are summarized in Ta- 
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Figure 4. {U - R), (B - R), and {V ~ R) colours of RS Oph 
versus R band magnitude for 2008 July 6. 



ble[l]and plotted in Figs[l][2]and[3l For each run we measure 
the minimum, maximum, and average brightness in the cor- 
responding band, plus the standard deviation of the run. 



3 COLOUR VARIATIONS 

As can be seen in Figs. [T] [2l [3] and Table 1, during our 
observations RS Oph exhibited variability on a time scale of 
1-30 minutes with amplitude 0.3 mag in V. The amplitude 
increases to shorter wavelengths and decreases to longer. 



3.1 Magnitude - Colour relation 

In Fig|4] we plot the colours of RS Oph versus the R band 
magnitude for 2008 July 6. The R band is chosen because we 
have the shortest exposure time. In order to relate magni- 
tudes in different bands, taken at non-matching start times 
and different exposure times, we linearly interpolated the 
light curves. The colours have been calculated over 3 grids 
with time resolutions 70, 120 and 360 seconds, to match the 
poorest filter sampling in the corresponding colour relation: 
70 s. for {V - R) vs. R, 120 s for {B - R) vs. R, and 360 s 
for {U — R) vs. R. Linear fits (of type y — a + bx) to the 
data points in FigUJgive: 



(U - R) = -9.73(±1.15) + 1.17(±0.11)i? 
{B -R) = -5.91(±0.71) + 0.79(±0.07)i? 
{V -R) = -0.86(±0.29) + 0.17(±0.03)i? 



(1) 
(2) 
(3) 
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Figure 5. Dereddened fluxes of the flickering light source of RS Oph. The solid line represents a black body flt. 

a) 2008 July 6 (circles), T^b = 9492 K, radius R = 2.6 Rq, located at distance d = 1.6 kpc. 

b) 2009 July 21 (diamonds) and 2009 July 23 (pluses) plotted together. The flt is Tjj, = 9370 K, radius R = i.T Rq. 

c) All data normalized to V band and plotted together. The flt is T;,;, = 9236 K. Symbols are identical to those in a) and b). 



The errors of the coefficients are given in brackets. These 
relations are obtained on the basis of the observations from 
2008 July 6. They are valid over the range 10.20 <R< 10.35 
mag. The Spearman's (rho) rank correlation gives p = 0.85 
(significance S.IO"*') for Eq.l, p = 0.66 (I.IO"'') for Eq.2, 
p = 0.57 (2.10~^°) for Eq.3. The significance in Eq.l-Eq.3 
is always << 0.001 indicating that all these correlations are 
highly significant and changes in the colours of RS Oph are 
correlated with brightness variations. 



3.2 Colours of the flickering source 

Bruch (1992) proposed that the light curve of CVs can 
be separated into two parts - constant light and variable 
(flickering) source. We assume that all the variability is 
due to flickering. In these suppositions the flickering light 
source is considered 100% modulated. Following these sug- 



gestions, we calculate the flux of the flickering light source as 
Ffi — Fav — Fmin , whcrc Fav IS the average flux during the 
run and Fmin is the minimum flux during the run (corrected 
for the typical error of the observations). Ffi has been cal- 
culated for each band, using the values given in Table 1 and 
Bessel (1979) calibration for the fluxes of a zero magnitude 
star. 

Adopting these results, we find that the flickering light 
source contributes about 6% of the light in the R band, 7% 
in V, 10% in B, and 13% in U. 

Following Snijders (1987) we assume interstellar extinc- 
tion Eb-v = 0.73 and an extinction law as given in Cardelli 
et al.(1989). Using the colour relations (Eq.l-Eq.3, derived 
in Sect. 3.1) in the interval R = 10.28 — 10.34 mag (in other 
words average 7? = 10.28 and we calculate average values for 
other bands, minimal brightness in 7? = 10.34, and we de- 
rive minimal fluxes in other bands), we obtain colours of the 
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flickering light source in RS Oph as {U-B)o = -0.67±0.06, 
{B - V)o = 0.08 ± 0.06, {V - R)o = 0.28 ± 0.08. On 2009 
July 21 we obtain (B - V)o = 0.22 and {V - R)o = 0.22. 
Using the entire runs in U band from 2009 July 23 and B 
and V from 2009 July 21 (as plotted in Fig[5]) we calculate 
{B - V)o = 0.22 and {U - B)o = -0.57. The values are 
similar to the colours of the flickering source in 1983 July- 
August {U~B)q = -0.71 and {B-V)o = 0.07-0.11 (Bruch 
1992). 



4 FLICKERING 

4.1 Temperature and size of the flickering source 

The derived {U — B)o colour corresponds to a B4V star 
(Te// ^ 17000 K) and a black body with T ^ 9000 K. The 
{B — V)o colour corresponds to an A5V {Teff ~ 9000 K) star 
and a black body with T ~ 10000 K. {V — R)o corresponds 
to a F8V star (T^ff = 7000 K). These estimates give an 
approximate temperature of the flickering light source Tfi — 
9000 - 12000 K. 

For the flickering light source we obtain magnitudes 
(corrected for interstellar extinction): 

2008 July 6 : U = 11.40, B = 12.06, V = 11.92, R = 
11.55 mag; 

2009 July 21 : B ^ 10.73, V = 10.52, R = 10.29, / = 
9.94 mag; 

2009 July 23: U = 10.00, V = 10.76, / = 10.11 - 10.45 
mag. 

The errors are in f/ ± 0.04, in B ± 0.05, in V ± 0.08, in 
i?±0.10 in /± 0.12. 

We adopt d = 1.6 ± 0.3 kpc as given in Bode (1987). In 
Fig[S]we plot these magnitudes transformed to fluxes. Using 
a black body fit {nfitld routine of IRAF), we calculate for 
the fiickering light source: Tji = 9492 ± 300 K; Rji = 2.6 ± 
0.3 Rq; Lf,r~.50 Lq (2008 July 6) and T;i = 9370 ±300 K; 
Rfi = 4.7 ±0.3 Rq; Lfir^ 150 L© (2009 July 21 and 23). In 
FigEJ; we plot all points normalized to the corresponding V 
band fiux. The fit gives Tfi = 9236 ± 200 K; 

Our /-band data from 2008 July 6 have considerably 
larger observational errors due to a technical problem which 
gave as a consequence a variable dark current pattern of the 
CCD. U from 2009 July 21 and B fi-om 2009 July 23 are 
short. They are therefore not used in the black body fits nor 
plotted in FigO 

There are several different estimates of the mass accre- 
tion rate in RS Oph. An estimate of Mace ~ 2.10"^^ Mq yr~^ 
arises from the required mass accreted on to the WD to 
give rise to the outburst plus the latest inter-outburst time 
scale of 21 years (Osborne et al. 2010). Assuming Rwd = 

0. 003 Rq (Starrfield et al. 1991), this results in a total ac- 
cretion luminosity Lace ~ GM^aM / R^d ~ 2900 Lq. Alter- 
natively, Walder et al. (2008) explore the hydrodynamics of 
wind accretion in the system and find Mace ~ 1O~*M0 yr~^, 

1. e. Lace ~ 146 Lq (which compares to Lace ~ 190 — 740 Lq 
from lUE observations at quiescence - Snijders 1987). The 
spectral energy distribution of the accretion disk around the 
WD is similar to mid-B spectral type (Dobrzycka et al. 1996) 
with luminosity 100 — 600 Lq . In practice therefore, the dif- 
ferent estimates give Lace ~ 140 — 2800 Lq, which means 
that 5%-35% of the total accretion luminosity is emitted by 
the flickering source (in other words Lfi/Lace ~ 0.05 — 0.35). 
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Figure 6. Position of the fiickering light source on a ((7 — B) vs 
{B — V) diagram. The dashed line is the main sequence. The solid 
line is a black body, (blue) triangles — CVs (from Bruch 1992), 
(black) squares - symbiotic stars (MWC 560, V407 Cyg), RS Oph 
is marked with (red) plus symbols, T CrB - (red) squares. 

4.2 Time delay 

We searched for time lags between the light curves in differ- 
ent bands, using the interpolation cross-correlation method 
(Gaskell & Sparke, 1986), but found no delays larger than 10 
sec, which is shorter than our time resolution (the exposure 
times). 

4.3 Difference between symbiotics and CVs? 

Fig.|S]represent a 2-colour diagram, {U — B)o versus {B — V)o 
for the flickering light source in a number of CVs and sym- 
biotics. The data for the CVs and T CrB are from Bruch 
(1992). We also plot 2 points for MWC 560 and 2 for 
V407 Cyg, which are calculated from the data of Gromadzki 
et al. (2006) 

The crosses represent RS Oph. The three crosses refer 
to the three epochs of observations - 1983, 2008, 2009. They 
lie well outside the area of the CVs. To address the possible 
errors of our measurement, we separated our light curves of 
RS Oph (plotted in Fig. [H into two parts (from UT 19:46 
to UT 20:26 and UT 20:27 - 21:22). We obtained for the 
flickering source {U - B)o = -0.61 {B - V)o = 0.11 for the 
flrst time interval, and {U - B)o = -0.74 {B - V)o = 0.15 
for the second. This indicates that for RS Oph the errors 
are approximately equal to the size of the symbols. 

To search for differences between the flickering of 
symbiotic stars and CVs we performed a two-dimensional 
Kolmogorov-Smirnov test (Peacock 1983; Fasano & Frances- 
chini 1987), using the data plotted on Fig|6] We compare 
colours of the flickering of the recurrent novae RS Oph and 
T CrB with those of the CVs. The test gives a probability 
of 2 X 10"'' that both distributions are extracted from the 
same parent population. 

We also performed the same test but comparing CVs 
with all symbiotic stars. The test gives a probability of 
4 X lO"** that both distributions are extracted from the same 
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parent population. The number of points is not high, but suf- 
ficient, as shown in Fasano & Franceschini (1987), and the 
difference is significant. This result indicates that in this di- 
agram there are clues to the cause of the differences between 
the flickering of these two classes of accreting white dwarfs. 

In the symbiotics of course, the mass donor is a red 
giant and the orbital periods are > 100 d. In the CVs the 
mass donors are late-type dwarfs and the orbital periods 
are ~ 1000 times shorter. To the best of our knowledge 
this is the first evidence that flickering differs in these types 
of accreting WDs and the physical cause of this difference 
warrants further investigation (see below). 



5 DISCUSSION 

In quiescence RS Oph varies irregulary between V — 10.0 — 
12.2 (CoUazzi et al. 2009, AAVSO light curves). This vari- 
ability is observed before as well as after the 2006 outburst. 
Darnley et al.(2008) detected an increase of the brightness 
in B from B = 13.5 to B = 11.9 mag and from V = 12.0 
to V — 10.5 for a year after the 2006 outburst. Worters 
et al. (2007) detected an increase in V from V = 11.3 to 
V = 11.9 for 2 weeks. The brightness of RS Oph at the 
time of our flickering observations is well inside these limits 
(unfortunately, there does not seem to be any published U 
light curve). The derived parameters of the flickering can be 
considered as typical for quiescence. 

After the 2006 outburst the flickering of RS Oph dis- 
appeared (Zamanov et al. 2006) probably as a result of (1) 
destruction of the accretion disk from the nova explosion or 
(2) a change in the inner disk associated with jet produc- 
tion (Sokoloski et al. 2008). The flickering resumed by day 
241 of the outburst (Worters et al. 2007). For the symbiotic 
star CH Cyg, Sokoloski & Kenyon (2003) observed changes 
in the flickering in association with the mass ejection event. 
Observations of CH Cyg and RS Oph therefore indicate that 
a connection does exist between the flickering behaviour and 
the ejection of matter from the WD. 

Different sites for the origin of the flickering have been 
discussed. They are all related to the accretion process: 

(i) the bright spot (the region of impact of the stream of 
transferred matter from the mass donor star on the accretion 
disk); 

(ii) the boundary layer (between the innermost accre- 
tion disk and the white dwarf surface); 

(iii) inside the accretion disk itself. 



(as in CVs), but via stellar wind accretion. As a result of 
the supersonic motion of the WD through the red giant wind 
there should be a shock cone and an accretion wake. Around 
the WD, it is likely that an accretion disk and/or a cocoon 
is formed. In the place where the matter flowing through 
the accretion wake encounters the accretion disk/cocoon a 
bright spot could be formed (analogous to the bright spot 
formed in the case of CVs, where the mass flow from the 
inner Lagrangian point encounters the accretion disk). 

(ii) Another possible site for the origin of the flickering 
in RS Oph is the boundary layer between the white dwarf 
and accretion disk. Bruch & Duschl (1993) estimated the size 
of the boundary layer (e) in RS Oph as e = 2.20. In their 
model Ld/Lu is connected with the size of the boundary 
layer {Ld is the the luminosity of the accretion disk, Lu is 
the luminosity of the boundary layer; see Fig 1 in Bruch 
& Duschl 1993). If we assume Ld + Lu ~ Lace ~ 110 — 
2000 Lq, and Lm ^ Lfi ~ 5 — 150 Lq, then we calculate 
Ld/Lu = 1.2—12, and the lower value agrees with the size of 
the boundary layer as estimated by Bruch & Duschl (1993). 
However, Tfi of RS Oph as derived in Sect. 4.1 is too low 
for a boundary layer. 

(iii) For UU Aqr, Baptista & Bortoletto (2008) found 
no evidence of flickering generated in regions (i) and (ii). 
They suggested that the flickering in UU Aqr is generated 
in the accretion disk itself and a possible reason can be tur- 
bulence generated after the collision of disk gas with the 
density-enhanced spiral wave in the accretion disk. 

The radial temperature profile of a steady-state accre- 
tion disk is: 



8naR^ 



1/2- 



(4) 



where a is the Stefan-Boltzmann constant, R is the radial 
distance from the WD. Using the parameters for RS Oph, a 
temperature Tfi ^^9500 K (the temperature of the flickering 
light source as derived in Sect. 14. l|) should be achieved at a 
distance R « 0.5 — 1 Rq from the WD. If (iii) is the place 
for the origin of the flickering of RS Oph, then it comes from 
R < 1 Rq from the WD. 

To understand more fully the nature of the flickering 
variability of RS Oph we need to acquire a set of spectra 
simultaneously with photometry and with time resolution 
~ 30 seconds (see also Sokoloski 2003). Such spectra poten- 
tially can directly give the spectrum of the flickering light 
source. 



The flndings of Bruch (2000) for HT Cas, V2051 Oph, 
IP Peg and UX UMa demonstrated that the flickering in 
these CVs can originate in both regions (i) and (ii). 

(i) The temperature and the size of the bright spot 
are derived for a number of CVs. A few examples of tem- 
peratures are: for OY Car Wood et al (1989) calculated 
black body T — 13800 ± 1300 K, and color temperature 
r = 9000 K; Marsh (1988) for IP Peg - T = 11200 K; 
Zhang & Robinson (1987) for U Gem - T = 11600 ± 500 K; 
Robinson et al. (1978) give T = 16000 K for the bright 
spot in WZ Sge. The temperature of the optical flickering 
source of RS Oph « 10000 K (see Sect. 4.1) is similar to 
the temperature of the bright spot for the CVs. The mass 
transfer in RS Oph may not be from Roche lobe overflow 



6 CONCLUSIONS 

We report our CCD observations of the flickering variability 
of the recurrent nova RS Oph, simultaneously with 4 tele- 
scopes in the U BVRI bands. RS Oph has a flickering source 
with {U - B)o = -0.62 ± 0.07, (B - V)o = 0.15 ± 0.10, 
{V - R)o = 0.25 ±0.05. 

For the flickering light source in RS Oph (1) we estimate 
Teff ~ 9500 ±500 K, which is similar to the temperature of 
the bright spot in CVs; (2) using a distance of d = 1.6 kpc, 
we find size Rji ~ 3.5 ± 0.5 Rq, and luminosity L/; ~ 
50 - 150 Lq. 

We find that on the {U — B) vs [B — V) diagram the 
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flickering of tlie symbiotic stars differs from tlie fiickering of 
the cataclysmic variables. 

The possible sites for the origin of the flickering are 
briefly discussed and it is suggested that time-resolved spec- 
troscopy should be carried out to explore this in more detail. 
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